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Abstract 

A general partial-wave multiple scattering theory for scattering from cylindrically symmetric 
potentials on a topological insulator (TI) surface is developed. As an application, the cross sections 
for a single scatterer and two scatterers are discussed. We find that the symmetry of differential 
cross section is reduced and the backscattering is allowed for massive Dirac fermions on gapped TI 
surface. Remarkably, a sharp resonance peak at the band edge of the gapped TI is found in the 
total cross section A tot , which may offer a useful way to determine the gap (as well as the effective 
mass of quasiparticles) on TI surface. We show that the interference effect is obvious in cross 
sections during the quasiparticle scattering between the scatterer pair, and additional resonance 
peaks are introduced in Atot when the higher partial waves are taken into account. 
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Topological insulator (TI) has attracted lots of interest in modern condensed matter 

fin nfi 

physics field since the advanced progress in theoretical [lH5] and experimental [6|-|9| studies of 
this new quantum matter phase. The ideal three-dimensional (3D) TI exhibits gapless modes 
on surface, and the backscattering is forbidden due to the time reversal invariance, which 
has been observed in the angle-resolved photo emission spectroscopy experiments on Bi 2 Te3 
surface 10] . However, if the Dirac fermions gain mass, a gap will be induced in system, 

llMl3|| and as a consequence, the scattering and 



which breaks the time reversal symmetry 



transport properties should be different from the gapless case. Due to the single Dirac- 
cone nature, no complicated intervalley scattering events occur on TI surface. Therefore, 
it can be expected that in situ measurement or even manipulation of impurity scattering 
processes on TI surfaces should pave a promising way to study electronic structures and find 



extraordinary quantum phenomena in various 
tunneling microscopy (STM) measurements 14 



I materials. Recently, a series of scanning 



16| and theoretical simulations 17|, |l8| have 



been performed on Bi2Te3 and Bii_ x Sb x surface states to study the scattering effect in the 
dilute impurity limit. On the other hand, when the impurities are located close to each other, 
the multiple scattering effects should be important. Especially, since the quasipartcle's spin 
is strongly coupled to its momentum, then the quantum interference between different spin 
states during multiple scattering process may display profound phenomena such as electric 
conductance weak (anti-)localization [l9| and Aharonov-Bohm effect 20( in STM signal. 
Therefore, it is clear that at present more deep multiple scattering studies are being called 
for. 

In this work we conduct a multiple scattering theory for the massive as well as massless 
Dirac fermions on surfaces of 3D TIs. A general partial-wave multiple scattering formula 
for scattering from short-range cylindrically symmetric potentials is developed. The cross 
sections for a single scatterer and two scatterers are discussed, and the interference effect 
during the quasiparticle scattering between the scatterer pair is found. We show that higher 
partial waves induce prominent corrections in cross sections, including additional resonance 
peaks. The symmetry of differential cross section is reduced and the backscattering can 
be observed for massive Dirac fermions due to the breaking of the time reversal symmetry. 
Especially, a sharp resonance peak at the band edge of the gapped Dirac fermion is found in 
total cross section, which may offer a useful way to determine the gap of the TI surface states. 
Note that our formulation is closely related to those given on semiconductor heterostructure 
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2l| and graphene [2 21-124} . Remarkably, unlike graphene, the TI surface is a single Dirac- 
cone system, which makes the present intravalley multiple scattering theory more natural 
and exact. 

We start from the effective-mass Hamiltonian near the Dirac point, 

H Q = hv F (er x k) ■ z + Acr 3 , (1) 

where -yF^5.0xl0 5 m/s is the Fermi velocity, er= (<7i, 02, 03) are Pauli matrices, k= (k x , k y ) 
is the planer momentum, and A=m*Vp is the band-gap, which is absent in the massless 
limit. The eigenstates of H are given by plane wave 

ik'f r jT } 

*o,±(r) = -^=(5 + , TiS-e i0 «) , (2) 



where 8±=yJ\e ± A|, and e= ± y (hvpk) 2 +A 2 with the upper (lower) sign referring 
to the electron (hole) part of the spectrum. The Berry phase of system given by 
— i J Q 27r d6f. (\&o,±| de k |^o,±) =^]7p 7r indicates that the backscattering is allowed (prohibited) 
if band-gap A^O (=0), which will be observed in differential cross section in the following 
discussion. To obtain the scattering theory from localized, cylindrically symmetric scatter- 
ers, it is convenient for boundary condition treatment to express the eigenstates for Eq. (CQ) 
in cylindrical coordinates, which are written as 



XS^ (r) = J±e«> ( MfP\ ±5-<f >e» f , 0) 

(1 2) 

for outgoing and incoming cylindrical waves about r=0, respectively. Here H( ' are Hankel 
functions of order / with variable kr. 

The potential of each symmetric scalar scatterer can be expressed as 

V (r) = V 6 (a - r) , (4) 

where a is the radius of scatterer. The incident plane wave centered about a single scatterer 
located at r n is given by 

n e^^V ( x g> (pj + X fX (Pj) , (5) 

i=— 00 

where X/± an d Xj ± represent outgoing and incoming waves about r n , with p n =r—r n and 
e tfln— Pni£tM _ By using the boundary condition at p n =«, we can get the scattered wave 



function 

oo 

^ (r) = s G +0 T -®™ + [siG+iTi + 8-iG- t T+] (6) 

1=1 

where the Zth-partial wave t— matrix is T ; ± =diag( P^, =B-P^i ) with 
P±=i(d x ±zd y y=^(d r ±i-d e ) 1 , and 



+l ~ 2|e| I i^^We^- £^) e ^ 
The variable of Hankel functions is kp n . The scattering amplitude of the Zth partial wave is 



= i+Jf (fcq) J m (fc'q) - i" Jj+i (fcq) J/ (A/o) 
A-if^ (fca) J z (fc'a) - A+H[ l) (ka) J l+1 {k'c N 



where A ± - 



k'—y (e— Vq) —A 2 /Hvf, and J; is the Bessel function of order /. In 



e-V .±A 

the massless limit, i.e., A=0, the above Eq. ([6]) can be reduced into a more compact form 
due to the symmetry of the band, which is written as 

n (r) = £ ^fi Gl (r, r n , e) T, , (9) 



k 

1=0 



where 



iW„ zrffW e -i(l+l)0n 

Gi oc I ' l+ r \ I , (10) 

and TJ=diag( p~ ; p+ ). Here, we have used the relation s_(j+i)=sj for A=0. 

It is easy to extend the above theory to multiple scattering problems of massive quasiparti- 
cles on TI surface. The total wave function for N scatterers located at positions r 1; r 2 , ■ • • 
is given by 'f (r) =$™+\l/j_ c , where 

\P± (r) = G (r) SM' 1 ^. (11) 

Here, G(r) is a 2x2N (2/ max +l) matrix, 

G(r)= (GCr.n), G(r,r 2 ), G(r,r N )) (12) 
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with G=[G +0 , G+i, G_i, ■ ■ • , G +Zmax , G_ imax ]. 5 is a 2iV(2/ max +l) x2iV (2Z max +l) diag- 
onal matrix with the nontrivial elements so and s±i. M is expressed by 



M 



GS, 



(13) 



where 



/ 



G 





G(2,l) 



G(l,2) 




G(1,N) 
G(2,N) 



\ 



\G(N,1) G(N,2) ••• J 



(14) 



with G(i,j) are 2 (2/ max +l) x2 (2Z max +l) matrices constructed by [Gf, 
(0</, l'<! max ). Finally, can be written as a 2iV (2/ max +l) xl vector, explicitly, 

T 



where 



(15) 



The 



T " [$ (r 4 )] , Tr [$ (r 4 )] , T+ [$ (r,)] , • • •, T" ax [$ ( P< )] , l£_ [$ fa)] 
multiple scattering problem will be simplified in the massless limit A=0, and similar to 
the case of a single scatterer, some transformations, including s_({ + i)=S{, (TJ + , 7] - )— )-T}, and 
(G+i, G-i)— 7-G;, should be performed in Eqs. flllH15p . Consequently, matrices M and 5 are 
reduced to 2N (/ max +l) x2iV (/ max +l) ones, which are similar to the problems of intravalley 



24J, so that we do not show the explicit 



multiple scattering of quasiparticles in graphene 
formulas for the case of A=0 herein. 

The above theory allows one to solve the multiple scattering problems on gapped or 
gapless TI surface with higher partial waves, which may be important as the distance between 
scatterers decreases as well as the scattering potential is strong. If the potential is very weak, 
s— wave (Z max =0) is enough. To obtain the scattering amplitude and the cross section, we 



apply the following approximation on Hankel functions in Eq. (fl2"|) for r— t-oo 



(k Pn ) _> 



(16) 



i 2l+1 7tkr 

and e ie "=p n - (x+iy) / 'p n —>f- (x+iy) =e lip , where f— (cos <p, sin <p) is a unit vector in the di 
rection of r. Finally, the scattered wave function can be approximated as 

1 / \ T e ikr 



^ (r) / (k, ip) — == ( 6+, TiS-e iv > 



(17) 
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FIG. 1: (Color online) The normalized differential (a) and the total (b-c) cross sections for a 
single scatterer. In (a), the red (blue) curve is for A=0 meV (A=50 meV), and e=320 meV. (b) 
is for A=0 meV, while (c) is for A=50 meV. Insert in (a) is a zoom of dA/dip in angle range 
[2ir/3, 47r/3], and the green dots indicate that the backscattering is allowed (forbidden) in the case 
of A 7^0 (A=0). In all calculations, Vq=2 eV, and a=l nm are chosen. 

where / (k, (p) is nothing else but the scattering amplitude. Therefore, the differential cross 
section should be 

dA/d<p=\f(k,(p)\ 2 , (18) 

and the total cross section is 

At* = ^ d V \f (k, <p)\ 2 = C A ^lm [e-^f (k, <p = 0)], (19) 

The second equality for A tot expresses the 2D optical theorem in terms of our definition of 
/. Ca= (<5++5 2 _) /2 |e| arises from the band gap, which reduces to unity in the massless 
limit. On the other hand, the scattered current j sc —^fVF ^nr- ^ fc ^ [cos <px+ sin tpy] , 
and the incident current 3i n = -f v F^j^- [cos 9kX+ sin 6ky], which result in 
dA/dip=r\j sc \ / \j in \ = \f (k,<p)\ 2 . In the following calculations, without losing gener- 
ality, we just consider the incident wave $™=-^=( ) T propagating along the 
positive x direction. 

Firstly, we discuss the single nonmagnetic scatterer problem. The numerical results of 
scattering cross section are shown in Fig. [TJ From Fig. Ufa), one can observe that the 
backscattering at potential is absent when the band gap decreases to zero, i.e., the differential 
cross section dA (A=0) /d<p\ w =0 [see the red curve in Fig. EJa)]. This is also can be seen 
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from the /th— partial wave scattering amplitude 



n 



A=0 




cos [(2/ + l)v?/2]e 



itp/2 



(20) 



which results in f L A (ip=n) =0. However, on the band-gapped TI surface, the Ith. partial 
wave scattering amplitude becomes 



try is broken by the nonzero gap in the surface Dirac spectrum. The blue curve in Fig. [T]^a) 
is for dA/dip corresponding to A=50 meV. It is clear that the symmetry of dA (A^O) /dip 
is reduced and the backscattering is allowed. This can be more clearly seen from the insert 
in[]Ta), which enlarges dA/dtp in angle range [2tt/3, 47r/3]. 

The total cross section A toi as a function of energy e is shown in Fig. \Hjo) arid [Q^c) 
for the gapless and gapped TI surface. On average, the contributions from higher partial 
waves are much smaller than s— wave (|/o| ^> \fi^o\) hi low-energy scattering regime, see the 
black curves in Fig. QJb) and[T]^c) which only takes into account the s— wave in calculations. 
However, we find that higher partial waves may introduce some remarkable corrections. For 
example, we find an additional resonance peak (cZ — like peak) around e=233 meV when the 
/j=2° partial wave is taken into account, while much higher partial wave has little corrections 
on A tot , see the red dashed curves in Fig. [TJ^b). This additional peak arises from the pole of 
| s/ = 2 1 , which could be obtained by expanding |s/=2| at this point. The pole of |s;=i| is out of 
the energy range we considered here, therefore, the corrections from p— wave is small in Fig. 
QXb). Furthermore, we can find that A tot for A^O [Fig. 0Jc)] is totally different from that 
for A=0. A sharp peak at the band edge e=A is found. In fact, e=A is a pole of for 
arbitrary Ith partial wave. Therefore, we can conclude that this sharp peak in A to t, which 
can be easily measured in scattering experiment, may provide a useful way to determine the 
band gap (as well as the effective mass m* of quasiparticle) of TI surface. The contributions 
from higher partial waves are obvious when e>A [the red dashed curve in Fig. [1(c)] . 

Now let us discuss the double-scatterer case. The analytical expressions of the scattering 
amplitude for multiple scatterers are complicated, especially when taking into account the 
higher partial waves and the band gap. In the limit of A=0, the differential cross section for 




(21) 



which indicates dA (A^O) /dtp\ = ^=7 f±t° (^ =7r ) ^° 

since the time-reversal symme- 
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FIG. 2: (Color online) The same as in Fig. [TJfor two identical scatterers separated by distance 
d=6 nm on a TI surface. r\={— 3, 0) and r2=(3, 0). Other parameters are chosen as the same as 
in Fig. [TJ 



two identical scatterers located at r\/2— (±<^/2, 0) becomes symmetric since M is symmetric. 
For instance, if we just consider the s— wave, we can obtain the scattering amplitude 

2 



A=0 





ink 



s 



E 1 

n,m=l 



2ra-l,2m 



2n,2m-l 



I 2n,2m 



(22) 



where 



/ 



Mo" 1 



1 

D 



(23) 



1 s H ±s #i \ 

1 TsoH! s H 

s H Ts H! 1 
\±s H 1 s H 1 

with D= v /det|Af | = l-sg{[^ 1) (fcd)] 2 +[#i (1) (^)] 2 }- It is easy to find that f (k, tt) =0, 
which indicates that the backscattering is also prohibited. This property can be clearly seen 
from the red curve in Fig. [2(a) as well as the inset of Fig. [2(a). However, the symmetry 
is reduced by the nonzero gap and the backscattering can occur, see the blue curve in Fig. 
12(a)- Besides, we find that the scattering along the directions of (/?~60°, 120°, 240°, and 
300° is forbidden since the interference effect along these directions is destructive during 
the multiple scattering process. In the forward direction (ip=0), however, the interference is 
still constructive. Moreover, comparing to the single-scatterer case, the width of resonance 
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FIG. 3: (Color online) Optical theorem Eq. (|19p for three scatterers. The black (red dashed) 
curves are results in the first (second) equality in Eq. (|19p . A=0 and / m ax=5 are chosen in 
calculations. 

peak is narrowed and the peak is pushed to lower energy due to the interference during 
quasiparticle scattering between the impurity pair. Also, the higher partial waves introduce 
additional resonance peaks in A tof , see the red curve in Fig. Mjo) and|2](c). 

In principle, one can get the cross section as well as other transport properties for 
N scatterers on the TI surface by the above theory. Before ending this paper, we calculated 
A tot for three scatterers to demonstrate the optical theorem in Eq. (fT9|) . For simplicity, we 
set A=0 and Z max = 3 in calculations, and the results are shown in Fig. [31 in which the black 
curve is obtained by numerical integral while the red dashed curve is obtained from the 
second equality in Eq. (fT9l) . We can see that they are in good agreement with each other. 
Therefore we believe the optical theorem in Eq. ( fT9l) should be proper to characterize the 
general multiple scattering process. 

In summary, we have developed a general partial-wave multiple scattering theory for 
scattering from cylindrically symmetric potentials on TI surfaces and applied to cross sec- 
tion calculations. We have found that higher partial waves may introduce some important 
corrections in low-energy scattering. Importantly, a resonance peak at the band edge of the 
gapped TI surface was found, which should be easily measured in scattering experiments. 
This prediction may provide a useful way to determine the band gap of TI surface states. Our 
formulation can be extended to magnetic impurity scattering as well as spin-orbit coupled 
scattering on TI surface. 
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